Schrodinger cat animated on a quantum computer 
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We present a quantum algorithm which allows to simulate chaos-assisted tunneling in deep 
semiclassical regime on existing quantum computers. This opens new possibilities for investigation 
of macroscopic quantum tunneling and realization of semiclassical Schrodinger cat oscillations. Our 
numerical studies determine the decoherence rate induced by noisy gates for these oscillations and 
propose a suitable parameter regime for their experimental implementation. 
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Since 1935 until recently, the metaphor of Schrodinger 
cat oscillations between life and death p| was considered 
as a purely theoretical concept. However, during the last 
decade such oscillations in a quantum limit were observed 
for two states of a Rydberg atom in a quantum cavity |^] 
and an experimental evidence was presented for a quan- 
tum superposition of macroscopically distinct states in a 
superconducting quantum interference device (SQUID) 
||. Manifestations of macroscopic quantum tunneling 
(MQT) were also observed in magnetization experiments 
with spin-ten molecular magnets Fe^ and Mriyi [^,^). 
In addition to their fundamental interest these experi- 
ments promise to provide important applications, e.g. for 
solid-state qubit realization || and information storage 
based on the Grover algorithm Q. 

The regime discussed in |lj assumes that the quan- 
tum tunneling takes place for a semiclassical object with 
a regular dynamics in two symmetric regions of phase 
space. Recently, the investigations of quantum chaos led 
to an extension of this concept to the phenomenon of 
chaos-assisted tunneling between islands of regular inte- 
grable dynamics separated by a chaotic sea ||. In this 
case, due to chaos the period T u of tunneling oscillations 
becomes very sensitive to the variation of system parame- 
ters and statistical description should be used to describe 
the average distribution of T u ||-[ll]] . This unusual tun- 
neling is strongly influenced by complex instanton orbits 
and scarring effects 10 and a chaos enhancement of 
tunneling rate by orders of magnitude may be reached 
by a small variation of parameters. The first direct ex- 
perimental observations of the chaos-assisted tunneling 
have been realized recently with cold fll2fl and ultracold 



atoms from a Bose-Einstein condensate |13 thus open- 
ing new possibilities for investigation of this interesting 
process. However, quantum tunneling is a very sensitive 
phenomenon and experimental studies are complicated 
by the decoherence produced by the environment. Due 
to that theoretical investigations of decoherence effects 
on MQT were initiated from the very beginning [ ^4|p5| 
and are continued up to now [ ^6p7| ] . 

In this Letter, we show that quantum entanglement 



and quantum computer simulations fll8|| can be efficiently 
used to study quantum tunneling in deep semiclassical 
regime. We illustrate this on the example of a quan- 
tum symplectic map (double well map) which has a rich 
phase space structure with integrable islands surrounded 
by chaotic sea. Our algorithm has certain similarities 
with the algorithms for the kicked rotator jl9| and for 
the sawtooth map pc| . It uses the quantum Fourier 
transform (QFT) ]2l| and simulates the dynamics of a 
quantum system with N levels in 0((log 2 A) 4 ) opera- 
tions per map iteration while any known classical algo- 
rithm requires at least 0(N log 2 N) operations. Only one 
work space qubit is required for computations so that n q 
qubits describe a physical system with TV = 2 rlq ~ 1 lev- 
els. Contrary to |2(]] the present algorithm simulates the 
quantum dynamics with mixed (chaotic/integrable) clas- 
sical phase space and is optimal for the investigation of 
chaos- assisted tunneling in semiclassical regime. Indeed, 
while for MQT in molecular magnets |4|j5|,|22]] the effec- 
tive Planck constant Ti is inversely proportional to the 
number of spins (Ti oc l/n q ) in our algorithm h cx 2~ rlq . 
Hence with only 10 qubits the algorithm allows to study 
MQT with the semiclassical parameter being larger by 
almost two orders of magnitude. At present the nuclear 
magnetic resonance (NMR) quantum computer can per- 
form QFT {H and operate with up to 7 qubits |J]. This 
allows to study the chaos-assisted quantum tunneling on 
existing NMR based quantum computers and to obtain 
important information about decoherence effects in MQT 
regime. 

In the classical limit the dynamics of our model is de- 
scribed by the double well map given by 

p = p — KdV(x)/dx , x = x+p (mod 2tt). (1) 

Here (p, x) are momentum and coordinate conjugated 
variables (— n < x < tt), the bars denote the variables 
after one map iteration and V(x) — (x 2 — a 2 ) 2 . In the 
limit K — > the map gives the one-dimensional inte- 
grable dynamics in the double well potential V(x) with 
the frequency of small oscillations ujq = 2\j2K. However, 
for K > the higher harmonics of finite step iterations 
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lead to appearance of chaotic component surrounding the 
stability islands located at x = ±a. A typical example of 
mixed phase space is shown in Fig.l 




FIG. 1. Poincare section for the double well map (|l|) at 
K = 0.04, a = 1.6, one chaotic and two regular orbits are 
shown in the cell (— 7r < x,p < n). 

The quantum evolution on one map iteration is de- 
scribed by a unitary operator U acting on the wave func- 
tion ip: 



ip = Uip 



-ip 2 /2h -iKV(x)/h 



(2) 



where p = —ihd/dx and ip(x + 2-K) = ip(x). In the follow- 
ing we take the dimensionless Ti — At: /N that corresponds 
to the case of quantum resonance |2q 1 with two classical 
cells (e.g. as in Fig. 1) on a quantum torus containing 
N levels. The semiclassical regime corresponds to Ti -C 1 
with discretized momentum p — Tin where n is an integer. 
The most efficient known classical algorithm simulating 
the quantum dynamics (g) is based on forward/backward 
fast Fourier transform (FFT) between p and x represen- 
tations. For a system with N levels it requires two FFT 
and two diagonal multiplications in p and x representa- 
tions and can be performed in 0{N log 2 (A r )) operations. 

The quantum algorithm simulates one map itera- 
tion (g) with N = levels in 0{nf) quan- 
tum gates operating on n q qubits with one qubit 
used as a work space. The initial wave function 
in x representation is coded in the physical register 
with n 

qubit. The action of kick Ut = ex-p(—iKV(x)/h) 
is diagonal in this representation and the simultane- 
ous multiplication of all N coefficients can be done 



q — 1 qubits in equidistant discrete points x m 
= X)m=o a ml m ) |0) with an empty work n q — th 



with a 3 \ 2 3 4 = or 1 . This step can be performed 
with approximately n q 4-qubit gates, namely control- 
control-control-phase shift (C (3) (/?)). The gate C^\p) 
is applied to each group of 4 qubits and transfers 
1 1111) to exp(-i/32^+-« + -' 3+ -' 4 )|llll) keeping other com- 
binations unchanged. Using the work qubit the gate 
C^ 3 ' (/?) can be expressed via two Toffoli gates T and one 
control-control-phase shift C^{(3) as Cf^^Jfi) = 

Tji ,j2,wC^j 3 j 4 ((3)Tj 1 ,j 2 ,w Here the indices indicate the 
qubits on which the gates do apply, f3 notes the rotation 
angle and w is the work qubit, which is reset to after 
C^ 2 '(/3). Thus, the action of kick is expressed via a se- 
quence of standard gates used for quantum computations 
|l]|6|. Indeed the Toffoli and C^{f3) gates can be ex- 
pressed via one and two qubit gates without addition of 
extra qubits Wq . The above computation is the most 




FIG. 2. (color) Time evolution of the Schrodinger cat: 
probability distribution W(x) at — tv < x < n is shown for 
different number of map iterations t, changing along y-axis 
from t = (top) to t = 180 (bottom). Here as in Fig.l 
K = 0.04, a = 1.6 and Ti = Aty/N with TV = 2 (n '- 1) . Quan- 
tum computation is done with n q — 6 qubits, ideal perfect 
gates (left) and noisy gates of strength e = 0.02 (right), and 
n g — 2090 gates per one map iteration (^) . At t = initial co- 
herent packet is located at x — —a. The color is proportional 
to the density: blue for zero and red for maximal density. 



in 3n q gate operations 
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then 
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difficult step in the algorithm and takes approximately 
2n q Toffoli gates and n q of C((3) gates for large n q . The 
multiplications by kick phases with lower powers of x 
are done in the similar way and require smaller num- 
ber of gates. After multiplication by Uk the algorithm 
exp(— i/3aj 1 aj 2 aj 3 aj 4 2 jl+j2 + j3+j4 ) is similar to the one used in |19 2(3]: the QFT changes 



Indeed, if x = Y^j=o a i^ : '> 
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x to p representation in 0(n 2 ) operations, the rotation 
Uh = exp(— ihn 2 /2) is realized in n 2 q of control-phase 
shift gates C(/3) and the backward QFT converts the 
wave function back to initial x representation. 

An example of the Schrodinger cat oscillations simu- 
lated by this quantum algorithm with ideal gates in the 
regime of chaos- assisted tunneling of Fig.l is shown in 
Fig. 2 (left). The time evolution of the probability distri- 
bution W{x), integrated over the work qubit, shows clear 
tunneling transitions between stability islands of Fig.l. 
The same evolution simulated by noisy gates is illustrated 
in Fig. 2 (right). Noisy gates are modeled by unitary ro- 
tations by an angle randomly fluctuating in the interval 
(—e/2, e/2) around ideal rotation angle. This noise intro- 
duces an effective decoherence rate F which destroys the 
tunneling oscillations after the time scale 1/T. 



the case of the coherent initial state. 
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FIG. 3. Probability for the Schrodinger cat to be alive W a 
(total probability for x < 0) as a function of time t for param- 
eters of Fig. 2. The time dependence allows to determine the 
period of chaos assisted tunneling oscillations T u — 90 and 
their decoherence decay rate F. Full curve shows the data 
without decoherece, points show the data for noisy gates with 
e — 0.01 (green) and e = 0.02 (black). The fit of data gives 
r = 1.9 x 10 -3 (dashed curve for e = 0.01). 



To determine the period of tunneling oscillations T u 
and their decay rate F it is convenient to analyze the 
time dependence of total probability W a (t) at x < (see 
Fig.3). The fit W a {t) - 1/2 oc e" rt cos(2nt/T u ) allows 
to obtain both T u and T. We note that the value of W a 
at given t can be obtained efficiently from few measure- 
ments that enables to determine T u and V. Moreover, the 
values of T u and V are not sensitive to the choice of initial 
0. As it was discussed in |28| for a similar 



state at t 

situation, this state should only have a sufficiently large 
overlap with the coherent state in the center of stability 
island. For example, the step distribution, W(x) = 2/N 
for x < and W(x) — for x > 0, which can be pre- 
pared efficiently, gives the same values of T u and V as in 
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FIG. 4. Dependence of decoherence rate of tunneling os- 
cillations r on the strength of gate noise e for different 
number of qubits n q — 6 (full triangles), 7 (diamonds), 8 
(squares) and 9 (circles). The selected map parameters are 
varied in the range 1.4 < a < 1.7, 0.04 < K < 0.06 at 
/2 (n "- 3 > that gave the tunneling period variation in 
the range 90 < T u < 1.8 X 10 5 . The straight line shows the 
average dependence (^). The data with noiseless work qubit 
at n q = 6 (K — 0.04, a = 1.6) are shown by open triangles. 
Logarithms are decimal. 

The dependence of the decoherence rate on the param- 
eters is shown in Fig. 4. The variation of T in a four orders 
of magnitude range is well described by the relation: 



r 



0.021e 2 n* 



(3) 



This relation can be understood using the following phys- 
ical arguments. For a given qubit, noise in each unitary 
gate gives a drop of the probability to be directed along 
the ideal direction by an amount of e 2 . Since at each map 
iteration the number of gates is n g ~ r& q the total decay 
rate is proportional to e 2 n g in agreement with Eq.([|). 
We note that similar estimates for the decoherence rate 
induced by noisy gates were also obtained for the Shor al- 
gorithm p9[| . At the same time the relation (3) is rather 
simple comparing to the decoherence rates discussed for 
MQT in SQUIDs @]l6) and molecular magnets @. 
One of the reasons for that is that the main step of the 
algorithm operates always with the same work qubit. If 
the unitary rotations of this qubit are noiseless then the 
decay rate T is significantly reduced (see Fig. 4). Hence, 
the quantum error correcting codes |l8| applied only to 
the work qubit can significantly reduce the decoherence 
rate, with a relatively small increase of the work space. 

Our algorithm allows to obtain interesting results 
about chaos-assisted tunneling even with a small num- 
ber of qubits. For example, the data of Fig. 3 can be 
obtained experimentally on the basis of techniques ap- 
plied in p3Lp4|]. The main obstacle for experimental im- 
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plementation of this algorithm is the decoherence. In- 
deed to observe the Schrodinger cat oscillations the de- 
coherence time scale 1/r should be much larger than 
the oscillation period T u . As it is usually the case for 
semiclassical tunneling, the later increases exponentially 
with the decrease of Ti. This implies very rapid growth 
of T u with n q : T u oc exp(S/h) oc exp(2™«5/87r), where 
S ~ 1 is a constant related to the classical action. For 
example, for K = 0.04, a = 1.6 the period changes from 
T u = 90 (n q = 6) to T u = 1.68 x 10 6 (n q = 9). There- 
fore even if the algorithm performs each map iteration 
(|J) in polynomial number of gates, exponential number 
of map iterations should be done to observe tunneling 
oscillations in deep semiclassical regime. Nevertheless, 
in the regime of chaos-assisted tunneling, the value of 
S can be easily varied H Q by changing the param- 
eters of the map {K and a) that allows to obtain not 
too large T u values for n q < 10, e.g. T u — 305 for 
K = 0.3, a = 0.5, n q = 10. Indeed, the value of S can be 
changed significantly by reducing the size of the stabil- 
ity islands embedded in the chaotic sea. In spite of the 
rapid growth of T u with the number of qubits n q the pro- 
posed algorithm uses them in an optimal way in order to 
reach the minimal effective Plack constant which scales 
as h oc 2~ 7lq . This situation is qualitatively different 
comparing to SQUIDs || where even for a macroscopic 
number of particles (analogous to n q ) the evolution is 
described by a Hamiltonian with two levels and effective 
% ~ 1 p^ ]. It also differs from the experiments [ p~2|J]~3| 
where effectively fi ~ 1 independently of number of cold 
atoms. 

In summary, our studies show that the Schrodinger cat 
can be animated in a deep semiclassical regime on exist- 
ing quantum computers |23j]24| ] with six or more qubits. 
Such experiments will give interesting information about 
the nontrivial regime of chaos-assisted tunneling in pres- 
ence of external decoherence. They will allow to deter- 
mine the effective accuracy of quantum computation, op- 
crability bounds and decoherence rates for the first gen- 
eration of quantum computers. 
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